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Abstract 

We present a self-consistent Schroedinger-Poisson scheme for simulation of electrostatic quantum 
dots defined in gated two-dimensional electron gas formed at n-AlGaAs/GaAs heterojunction. 
The computational method is applied to a quantitative description of transport properties studied 
experimentally by Elzermann et al. [Appl. Phys. Lett. 84, 4617 (2004)]. The three-dimensional 
model describes the electrostatics of the entire device with a quantum dot that changes shape 
and floats inside a gated region when the applied voltages are varied. Our approach accounts 
for the metal electrodes of arbitrary geometry and configuration, includes magnetic field applied 
perpendicular to the growth direction, electron-electron correlation in the confined electron system 
and its interaction with the electron reservoir surrounding the quantum dot. We calculate the 
electric field, the space charge distribution as well as energies and wave functions of confined 
electrons to describe opening of two transport channels between the reservoir and the confined 
charge puddle. We determine the voltages for charging the dot with up to 4 electrons. The results 
are in a qualitative and quantitative agreement with the experimental data. 

PACS numbers: 73. 21. La 



1 



I. INTRODUCTION 



A quest for a nanodevice that would store a quantum bit in the spin 1 ' 2 of a confined 
electron and allow for its manipulation is a main factor that motivates the research on gated 
electrostatic quantum dots. Realization of a spin quantum gate requires application of mul- 
tiple quantum dots coupled in way that can be controlled during the device operation. First 
electrostatic quantum dots were formed in a vertical configuration^^^ of a semiconductor 
heterostructure containing a single or multiple quantum well surrounded by a single gate 
creating the lateral confinement. The strength of the tunnel coupling between the vertical 
dot and electron reservoirs depends on the applied barrier thickness and composition, and 
therefore it is fixed for each device. Similarly, for vertical artificial molecules the interdot 
coupling is defined at the production stage. 

A full control of the interdot coupling is possible in quantum dots formed in gated two- 
dimensional electron gas.-^ 1 ^^^ 1 ^ These structures are produced by deposition of mul- 
tiple gates on top of the n-AlGaAs/GaAs structure containing a two-dimensional electron 
gas at the heterojunction. The system of gates is designed to locally deplete the electron gas 
in order to tailor a quantum dot of the surrounding electron reservoir. The voltages applied 
to the multiple gates define the confinement potential for the trapped electrons and control 
the tunnel barrier between the dot and the reservoir. Thus, the coupling of the confined 
artificial atom to the environment can be intentionally tuned from an open to a closed dot 
regime. For nearly open dots the Kondo effect and co-tunneling phenomena are observed.— 
Coupled dots with voltage tunable interdot barriers are also produced.— Quantum dots 
formed in the gated two-dimensional electron gas are used for investigation of spin depen- 
dent transport 1 — and confined spin relaxation.— A capacitive coupling between the dot and a 
quantum point contact defined in the same structure^&iiii 2 -^ is used to probe the confined 
states by the conductance measurements. The purpose of the present paper is to provide a 
theoretical description of a quantum dot formed in the gated two-dimensional electron gas. 
We focus our attention on the nanodevice that is probed by the quantum point contact as 
described in Refs. f^] and p^]. 

The electrostatic confinement potential in vertical quantum dots was thoroughly stud- 
ied in Refs. p&2Ii22j22i2ii2£i2iii21] . Less attention was paid to dots based on the gated two- 
dimensional electron gas. In particular, a theory for a double planar dotr^ was provided 
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in Refs. p£] and p 9 .]. Ref. p£] describes a structure of a triple quantum dot. Theoretical 
modeling of planar dots is for several reasons more difficult than modeling of the vertical 
structures. In vertical dots the electrons are confined inside a relatively deep quantum well 
with the lateral confinement strength controlled by a single gate. In planar structures the 
confinement is entirely due to the voltages that are applied to multiple gates and the formed 
potential cavity is typically shallow. Therefore, both the confinement potential and the 
few-electron eigenproblem have to be calculated with a special care. Moreover, as we show 
below, the gates not only fix the strength of the dot confinement but also change the shape 
of confined charge island which floats with the voltages inside the gated area. Interaction 
of the confined electrons with the reservoirs is also more complex due to variation of the 
depleted gas region with the applied potentials. The vertical dot^2^^^^^2^ often have 
circular symmetry,^ 2 which is not the case for the planar structure discussed below. 

II. THEORY 

A. Model of the structure 

The nanodevice^ii which we aim to describe is constructed on a basis of a planar semi- 
conductor heterostructure of GaAs/n-Alo.3Gao.7As in which the two-dimensional electron 
gas is created at the GaAs side of the junction. A cross section of the layer structure is 
presented in Fig. 1. On the substrate side there is a 1500 nm thick layer of undoped GaAs 
with a blocking AlGaAs barrier deposited on top. Lower part of the barrier (20 nm thick) is 
undoped and the upper (65 nm) is heavily doped with donors. On top of the barrier there 
is a 5 nm thin layer of n-doped GaAs. The donor states in AlGaAs stay 200 meV above 
the conduction band minimum of GaAs. Therefore, the electrons pass to the GaAs layer 
but stay localized under the barrier due to the Coulomb attraction by the charge of ionized 
donors. 

In our model we assume the shape of gates deposited on top of the structure (see Fig. 2) 
according to Refs. p&ik]. A negative voltage applied to the electrodes depletes the electron 
gas underneath and forms a lateral confinement potential in the center of the gated area. 
For properly adjusted voltages a few electrons stay localized in the center of the structure 
forming an artificial atom. The electrons are confined in the vertical direction by the barrier 
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FIG. 1: (color online) Structure of layers used for formation of a gated quantum dot in the two- 
dimensional electron gas (according to Ref. f^]). 

formed due to the GaAs/AlGaAs conduction band offset Uj,. The electrostatic potential 
0eist( r ) is responsible for the lateral confinement as well as for the potential that closes the 
dot from the substrate side. The quantum dot confinement potential is therefore expressed 
by 

^conf(r) = U h (z) - |e|&bt(r). (1) 

Ub equals zero in GaAs layer and the conduction band offset in the AlGaAs barrier. The 
electrons are additionally confined by an in-plane magnetic field of 10 T applied parallel to 
the surface of the layer structure.- 1 ^ Application of a strong magnetic field within the plane of 
confinement^ was previously discussed in context of tuning (reduction) of the electron tunnel 
coupling for vertical artificial molecules. To the best of our knowledge the Schroedinger- 
Poisson problem for the two dimensional electron gas with an in-plane orientation of the 
magnetic field was never solved previously. 

B. Sources of electrostatic potential 

The total electrostatic potential t ot( r ) is influenced by voltages applied between the 
substrate and metal electrodes on top of the structure as well as by the charge distribution 
inside the device. The charge density ptot( r ) is a sum of three contributions that have 
different nature and distribution 

Ptot(r) = p eq d(r) + Pd(r) + p e i(r). (2) 

The first contribution p eq d( r ) is the distribution of the charge confined in the quantum dot 
that is found by the solution of a few-electron quantum eigenproblem that is in the present 
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work obtained by the configuration interaction method. The second contribution Pd( r ) is 
the ionized donor space charge in the AlGaAs barrier, and the third p e i(r) is the charge 
density of the el — 
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FIG. 2: Geometry, size and position of the metal gates deposited on top of the semiconductor 
heterostructure (according to Ref. f^]). 

According to the superposition principle the total electrostatic potential ^tot(r) can be 
expressed as a sum of contributions of all the charge densities. We separate the potential due 
to the confined electrons eq d of the total potential. In this way we obtain a component of the 
potential <fi e \ st that enters the formula (jTJ) for the confinement potential of the electrostatic 
dot 

0totO) = 0eqd(r) + 0clst(r). (3) 

Potential of the dot-confined charge is calculated directly from the Coulomb law 

eqd (r) = j^-f P f^dv>. (4) 
47ree J |r — r| 

Potential e ist of all the other sources is found from the Poisson equation 

V 2 elst (r) = (5) 

The proper boundary conditions are naturally given for the total potential tot (r). The 
boundary conditions for the Poisson equation (jHJ) are obtained from boundary conditions 
that are known for the total potential 0tot(r) by a simple subtraction 

</W0) = 0tot(r) - </> C qd(r)- (6) 
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Charge density p(r) in Eq. (J5J) is a sum of the charge densities of ionized donors and of the 
electron gas 

P(r) = Pd(r) + Pci(r). (7) 

The donors are ionized in the spatial region where the total electron potential energy cal- 
culated with respect to GaAs layer is larger than the donor binding energy. For the donor 
energy level taken as the reference energy we obtain the ionization condition 

, x J for -M0 tot (r) < E D 

Pd(r) = < , (8) 

y |e|n D (r) for -|e|0 tot (r) > E D 

where nj)(r) is the density of donor impurities, and E D is the donor binding energy. As- 
sumption of a homogenous (continuous) donor distribution is justified by the presence of an 
undoped AlGaAs buffer (see Fig. 1). 

Calculation of p e i( r ) _ the second charge density entering Eq. fl7|) that corresponds to the 
electron gas confined at the hetero junction - is a nontrivial task. It can be exactly evaluated 
only in the asymptotic region, i.e. at a large distance of the electrodes. In the neighborhood 
of the electrodes one has to introduce an approximate treatment (see II. D). 

C. Potential and charge distribution in the asymptotic region 

At a large distance of the electrodes the electric field is parallel to the growth direction (z) 
and the electrostatic potential does not depend on the other two coordinates. The potential 
distribution results from an equilibrium between the ionized donors in the AlGaAs layer and 
the electron gas confined below. We choose the y axis as parallel to the external magnetic 
field. We adapt the Landau gauge 

A(r) = (-B*,0,0) (9) 

which leads to the Hamiltonian of an electron confined at the hetero junction 

_ cP_ cP_\ z d_ 

2m \dx 2 dy 2 dz 2 J c dx 

+^uj 2 c z 2 + U b {z) - \e\(f> e]Bt (z), (10) 

where lo c = and m is the electron band mass. Since this Hamiltonian commutes with 
momentum in both x and y directions its eigenfunctions are expected to be of form 

$ as (r) = Cexp(iqx)exp(iky)(j) as (z). (11) 



Eigenequation for ^(z) is obtained by substitution of Eq. (TTTj) into tTTUj) 

H{z)<t>™{z)=e nq <t>™{z), (12) 



where 



= "2^^ + _ Zo)2 + Ub{z) " i e i^W' ( 13 ) 



with zq = %q/muj c . 

Eigenvalues e nq and the eigenfunctions (j>nq( z ) are labeled by a natural quantum number 
n and depend on the wave vector q in the direction perpendicular to the magnetic field 
direction. Wave vector q enters the Hamiltonian operator (13) through shifted harmonic 
oscillator minimum Zq. The total electron energy eigenvalues are given by 

h 2 k 2 

E n kq = £nq + ~ • (14) 

2m 

The electrons confined at the hetero junction have energies below the Fermi energy. In the 
present calculations the Fermi energy is taken as the reference energy level. Only the states 
with E n kq < can be confined at the heterojunction. Given the H(z) eigenfunctions one 
calculates the charge density of the electron gas 

f%(z) = -2\e\ \^nkq(x,y,z)\ 2 (15) 

nkq 



^ I dkdq\^ q (z) 



^, , dq f F dk\^ q {z)\\ 

_ Z7r Js nq <0 J-k F 



where kp = \l — 2m f i nq . By integration over k we obtain 



PeTW = / W-^ICWI 2 - (16) 

„ 71 Je nq <0 v n 

Apart from the electron gas another source of the electric field is the ionized donor dis- 
tribution p d calculated according to Eq. (jSj) in which we identify = (p^ st since at an 
asymptotically large distance of the quantum dot potential </> eq d vanishes. Finally, the elec- 
trostatic potential is calculated from a single- dimensional Poisson equation 

92 .as , z) = Pgf(*)+Pd(^ 



dz^ elsty ~' ee 



Since p^f appearing in Eq. (17) is calculated with the Schroedinger equation (fT2l) . which 
in turn contains the electrostatic potential, both the equations are solved in an iteration 
till self-consistency is reached. Eq. ( I12p is solved below the barrier in a computational 
box long enough to allow the electron charge density ffTB]) to vanish before its end. The 
equilibrium solution is obtained when two conditions are met. The first one results from the 
charge neutrality which requires that the number of ionized donors is equal to the number of 
electrons trapped at the interface, which results in vanishing electric field at both ends of the 
region where the Poisson equation (ITT)) is solved. In fact, whenever a self-consistency of Eq. 
(12) and (17) is obtained the potential derivatives at both ends of the computational box 
vanish. The second equilibrium condition requires that the number of ionized donors is such 
that the variation of the electric potential on the entire space charge (on both sides of the 
hetero junction) equalizes the jump in the conduction band at the GaAs/AlGaAs interface. 
The latter is due to the fact that the electrons occupy all the states below the Fermi energy. 
This includes both the donor states in the barrier (not all the donors are ionized) and the 
states trapped below the barrier. 

The surface density of the electron charge accumulated below the barrier is obtained by 
integration of p^{z) along the growth direction 

/oo 
dzpf^z), (18) 
"OO 

with the charge neutrality condition 

d = » (19) 

where cxff is the surface density of ionized donors. The surface densities depend on the 
conduction band offset. Under the assumption of a homogenous donor distribution the 
asymptotic ionization range d fulfills the condition = n^d. The second equilibrium 
condition is obtained for a properly chosen d. For a nominal composition of the barrier 
Al x Gai_ x As with x = 0.27 the barrier height is Ub = 229 meV. For that value one obtains 
the asymptotic surface density of cr as = 3.5 x 10~ n cm -2 , which is close to the nominal 
experimental valued of 2.9 x 10~ n cm -2 . The deviation of the calculated density off the 
nominal value may result from the neglect of the exchange interaction in the electron gas. 
We decided to reduce the barrier height to Ub = 200 meV for which the calculated density 
is equal to its nominal value. 

8 




z [nm] 



FIG. 3: Electron potential energy (solid line) and the charge density of the electron gas (dashed 
curve) at the AlGaAs/GaAs heterojunction (GaAs is at the positive side of z). Thin horizontal 
dashed line shows the Fermi energy pinned at the donor impurity level in AlGaAs. E D is the donor 
binding energy that enters formula (8). 

The potential and the electron density calculated for the asymptotic region according to 
the procedure explained above is presented in Fig. 3. From the potential dependence on z 
we can see that the charge neutrality condition (zero electric field at both ends of the box) is 
fulfilled. We also notice that the electrostatic potentials of both sides of the junction differ 
exactly by the donor binding energy. This is because in the discussed structure the donor 
impurity level in the heavily doped AlGaAs layer defines the Fermi energy. 



D. Potential and charge distribution near the quantum dot 

Under the electrodes the electron potential energy is still positive also below the AlGaAs 
barrier. This removes electron gas from the region below the electrodes. The electron 
potential energy decreases with the growing distance of the gates and it eventually becomes 
negative in the region where the electron gas is not completely depleted. The spatial variation 
of the electron density is crucial for the shape and width of the potential barriers which 
separate the quantum dot of the electron reservoirs. An account for the electron dependence 
on the potential energy is taken in a following approximate manner. We assume that the 
electron gas density is zero wherever the local potential energy exceeds the Fermi energy 
(Ep = 0). In region where the local potential energy is negative we assume that the electron 
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density is proportional to its absolute value 



for r < 

p e i(x,y,z) = { , (20) 

rp™(z) for r > 



for 

-\e\(j) tot (x,y,z c ) + U b (z c 



-\e\^ ot (z c ) + U b (z c ) ' (21) 
with t < 1, where 0j^ t = 0gf st and z c is the center of mass of the asymptotic electron density 



-L J dzp™(z)z. (22) 



The adopted formula (20) simulates the depletion of the electron gas in the region of positive 
potential energy. In the asymptotic region r tends to 1, which guarantees that the known 
value of the electron density is found far away from the gates. 



E. Boundary conditions for Poisson equation in three dimensions 

Poisson equation (J5j) is solved in a three dimensional rectangular region which contains the 
quantum dot and a sufficiently large part of the electrodes. A standard test for the acceptable 
size of the box consists in performing the calculation in function of the box dimensions. We 
find that the results eventually saturate for a rectangular box of side lengths L x = L y = 600 
nm and L z = 400 nm. The position of the adopted box with respect to the electrodes is 
presented in Fig. 2. In the growth direction (z) the box covers 200 nm on both sides of the 
AlGaAs/GaAs heterojunction. 

A unique solution of the Poisson equation is obtained for boundary condition given on the 
surface of the computational box and on the metal electrodes that are inside the integration 
region. The boundary conditions are naturally given for <p tot . We calculate the conditions for 
0eist according to Eq. (jSj). On the electrodes the total potential is constant and determined by 
applied voltages. On the surface of the electrodes we assume Dirichlet boundary conditions 

0tot = U X + U s , (23) 

where X = P,T, M and R enumerates the electrodes (see Fig. 2), Ux is the applied voltage 
and Us is the Schottky potential which at the metal/GaAs contact is Us = —0.65 V. On 
the lateral sides of the computational box we apply Neumann boundary conditions 

n ■ V&ot = 0, (24) 
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where n is the vector normal to the surface. The Neumann conditions are consistent with 
the Gauss law for charge neutrality of the computational box content. On the vertical walls 
of the box (parallel to the growth direction z) this boundary condition is equivalent to the 
assumption that the electric field is perpendicular to the heterojunction, which agrees with 
the boundary condition used in the asymptotic region. 



F. Electrons confined in the quantum dot 

The potential minimum that is found in the central region between the electrodes traps 
several electrons provided that the applied voltages are not too negative. Hamiltonian for 
the system of N electrons with Landau gauge writes 

V ft*, .. . d 



( h 

i=l v 

in \ 
+y^c(^ - z o) 2 + U b (z) - |e|0eist(ri)J 



N N 2 

The eigenproblem Hn^n = En^n is solved with a configuration interaction approach in the 
basis of Slater determinants built of single-electron wave functions that are calculated with 
a finite-difference technique on a three-dimensional mesh. The convergence for the ground 
state energy of 4 electrons is achieved for the basis containing all the Slater determinants 
that can be constructed of the wave functions of 10 lowest-energy single-electron levels. 

The parameter z Q introduced by the gauge transformation is taken to minimize the total 
energy of N electrons. In our calculations we adopt z ~ 12 nm which coincides with the 
center of the density of the dot-confined electrons^ (see also Fig. 2). 

The few-electron wave function is used to evaluate the confined charge density 

Peqd(ri) = -|e| Jdr 2 dr 3 ...dr N \y N (r 1 ,r 2 ,r 3 ,...,r N )\ 2 . (26) 

The ground-state energies for iV and N — 1 electrons determine the electrochemical po- 
tential of the iV-electron quantum dot 

/ijv = E N — En-i- (27) 

The dot is filled with exactly N electrons when /xjv < Ep < //jv+i- Charging lines that are 
detected in the experiment correspond to ^ — E F . 
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G. Numerical procedure and self-consistency 

The Poisson equation (jSJ) is solved on a three dimensional mesh with a finite difference 
method. The adopted mesh steps Ax = Ay = 12.5 nm and Az = 2 nm are sufficient to 
describe the charge distributions and the shapes of electrodes. Smaller step in the growth 
direction is necessary because of a strong localization of the electron gas at the interface. 
Same step sizes are applied in the Schroedinger equation. 

According to expression (7) the charge density that enters the Poisson equation (jSJ) is a 
sum of charge densities of ionized donors that depend on the potential in a manner defined 
by Eq. (8) and the electron gas accumulated at the GaAs/AlGaAs interface. The latter 
depends on the total potential according to Eq. (20). Electrons trapped by the quantum dot 
are the third charge density [Eq.(26)] present in the nanodevice. The Schroedinger-Poisson 
calculations are iterated till the self-consistency of the three charge distributions with the 
total potential is obtained. The iteration accounting for the dependence of the electrons 
and ionized donors on the electrostatic potential requires application of an under-relaxation 
technique to ensure stability. The under-relaxed iteration is usually slowly convergent and 
requires several hundred thousands iteration for the entire mesh. The convergence is radically 
improved for a strategy of simulated cooling of the system. The measurements are performed 
at temperatures of several mK, for which occupation of electron states in Eqs. (8) and (20) 
has a nearly binary distribution. For non-zero temperatures we replace the formulas for the 
distribution of ionized donors and the electron gas used in the above theory by expressions 
accounting for the Fermi statistics 



which tends to distribution given by (8) in T = limit. For the electron gas at the interface 
we use the formula 



which in the limit T —>■ tends to Eq. (20) with r given in (21). We start the iteration at 
T = 15 K, in which the convergence is quickly reached. Next the temperature is gradually 
reduced to the nominal value. With the simulated cooling the convergence is obtained with 
an overall number of iterations that is reduced ten times. 




(28) 




(29) 
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III. RESULTS 



The numerical procedure described above was applied for the device described in Ref . f^] . 
Basic parameters of the nanodevice including composition and width of the semiconductor 
layers, the position and shape of the gates, and the applied voltages are adopted accord- 
ing to the experimental data.— The theory contains a single free parameter: the dopant 
concentration in the AlGaAs barrier, which is not known with a sufficient precision since 
the position of the transmission lines turns out to be extremely sensitive to its even small 
variation. The donor density appearing in formulas (8) and (28) was set to reproduce the 
charging of the dot with the first electron when the voltages are set as in the experiment 
(see below). 

The two-dimensional electron gas in the asymptotic region acts as an electron reservoir 
whose electrochemical potential is set by the voltage of the source and drain that are con- 
nected to the electron gas in the asymptotic region. We adopt the potential applied to the 
source, drain and the electron gas as the reference value for the voltages V$ — Vp — 0. 
Following the experiment^ we assume that voltages applied to the electrodes are: Vp = 0, 
Vp = 1.5 V, Vm = —1.07 V and Vp = —0.096 V. Under these voltages a first electron 
occupies the dot.— Figs. 4(c) and 4(b) present the electrostatic potential distribution in 
x,y plane for z = z a = 12 nm [used in Eq. (25)] calculated for np = 20 x 10 16 cm -3 and 
rip = 25 x 10 16 cm -3 , respectively. In both figures the blue line shows the zero of the electro- 
static potential. Note, that in Fig. 4(c) the zero potential is found far from the center of the 
device which indicates that the dot cannot trap any electrons. In Fig. 4(b) the zero level 
encircles quite a large region which turns out to trap several electrons and not a single one. 
We found that a single electron occupies the dot with the zero binding energy for the donor 
concentration np = 21.641 x 10 16 cm~ 3 . The corresponding potential profile in presented 
on a surface plot in Fig. 4(a) and as a contour plot in Fig. 5(a). The shaded region in 
Fig. 5(a) shows the electrodes. The potential has a well developed minimum of negative 
potential that is just enough to trap a single electron. The electron reservoir which encircles 
the dot is presented in Fig. 6 which shows a density map of the electron gas confined at the 
heterojunction. For the applied voltages the oscillating plunger gate (P) voltage takes the 
electron out of the dot to the reservoir and attracts it back with the electron tunneling along 
the transmission channel that is opened parallel to the x direction. The opened channel is 
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FIG. 4: Electron potential energy U(x,y) = — \e\<p e \ B t(x,y, z Q ) calculated at a distance of z Q = 12 
nm of the AlGaAs barrier for voltages V P = 0, V T = 1.5 V, V M = -1-07 V and Vr = -0.096 V. 
Assumed donor impurity concentration is no = 21.641 x 10 16 cm -3 in (a), rip = 25 x 10 16 cm -3 
in (b) and no = 20 x 10 16 cm -3 in (c). Blue line shows the U = contour. 

visible in Fig. 5(a) for x, y ~ 200 nm. In Fig. 5(b) and (c) the potential cavity is surrounded 
by potential barriers that are too deep for the electron to cross. 

Formation of a single-electron artificial atom is obtained also for other voltages. Lowered 
Vm voltage results in an increase of the electron potential energy in the center of the nan- 
odevice that can be compensated by an increase of Vr. A transmission line can be traced on 
the (Vm, Vr) plane. Similar transmission lines are observed for the tunneling of the second 
electron from and into the dot. In Fig. 7 we present the transmission lines calculated for 
1-4 electrons which are compared to the experimental results— We note that the calcula- 
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FIG. 5: (Color online) Equipotential contours and regions of negative electron potential energy 
(plotted in blue) for nn = 21.641 x 10 16 cm -3 at 12 nm of the barrier. In (a) the voltages are same 
as in Fig. 4. In (b) and (c) (-1.12V,-0.93V) and (-1.18V,-0.90V), respectively. Shaded regions show 
the metal gates. 

tions very well reproduce both the slope and the positions of the transmission lines. The 
experimental lines for the first and the second electron vanish for larger Vm, which is due to 
the increase of the potential barrier closing the transmission channel. For third and fourth 
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FIG. 6: Contour plot of the density of the electron gas at 12 nm of the barrier for the voltages 
used in Fig. 6(a), the darker shade of grey the larger density. 
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FIG. 7: Transmission lines as calculated (solid lines) and measured^ (thick dashed lines). Above 
the highest line the dot is empty. 

electron the lines reappear for larger Vm, which is explained by opening of a second trans- 
mission channel. Increase of the potential barrier and formation of the second transmission 
channel can be observed in Figs. 5(b) and 5(c) which present the potential distribution with 
a single dot-confined electron calculated for Vm,Vr voltages equal to (-1.12V,-0.93V) and 
(-1.18V,-0.90V), respectively. Fig. 5(b) shows the case where the left channel of Fig. 5(a) 
is closed. In Fig. 5(c) we see lowering of the potential barrier at right of the dot center. 
The tunneling probability is still not large enough to allow for an experimental observation 
of the transmission. 
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Interaction of the confined electrons with the electron gas leads to an appearance of 
nonlinear effects in the Schroedinger equation2&2i f or th e dot-confined electrons. The non- 
linear effects increase the electron binding energies and result in self-focusing mechanism for 
a single confined electron. An account of this self-focusing potential is taken in the Poisson 
equation (151) . The Poisson equation does not account for the energy relaxation of the dot- 
confined-electron system to the reservoir.— The energy relaxation leads to a decoherence of 
the dot-confined quantum states. The energy relaxation rate is likely to be significant due 
to small density of the electron gas in the neighborhood of the dot.— 



IV. SUMMARY 



We presented a theory describing phenomena appearing in a multielectrode device of 
gated two-dimensional electron gas containing a quantum dot. In the asymptotic region of a 
large distance from the gated area the calculation consists in solution of a single- dimensional 
Poisson-Schroedinger problem for the electron gas in a strong magnetic field parallel to the 
semiconductor surface. Solution in the asymptotic region is used to determine the electron 
gas density in the three-dimensional Poisson-Schroedinger problem for the gated region 
with a quantum dot. The energies and charge densities of several confined electrons were 
calculated with a configuration interaction approach. The presented theory includes a single 
fitting parameter - the donor concentration. All the other device parameters: the layer 
structure, the shape, size and position of the electrodes as well as the applied magnetic field 
vector are taken of the experiment jy&II The slope and positions of the transmission lines 
were reproduced with a very good quantitative agreement with the experiment. We discussed 
the electric field distribution in the device for voltages corresponding to transmission lines 
observed in Ref. [9] with a particular attention to the inhomogeneities creating the quantum 
dot confinement and the surrounding potential barriers that separate the artificial atom of 
the electron gas. 

For voltages corresponding to the transmission lines the weakest bound electron is stim- 
ulated to oscillate between the dot and the reservoir by the oscillating plunger voltage. We 
demonstrated opening and closing of two transmission channels in the barriers that allow 
for the oscillations of the confined charge. 
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